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Abstract 

We consider a phenomenologically viable 50(10) grand unification model of the 
unification scale Mg around 10^^ GeV which reproduces the MSSM at low energy and 
allows perturbative calculations up to the Planck scale Mp or the string scale Mgt- 
Both requirements strongly restrict a choice of Higgs representations in a model. We 
propose a simple 50(10) model with a set of Higgs representations {2 x 10 + 16 + 
16 + 45} and show its phenomenological viability. This model can indeed reproduce 
the low-energy experimental data relating the charged fermion masses and mixings. 
Neutrino oscillation data can be consistently incorporated in the model, leading to 
the right-handed neutrino mass scale Mr ~ Mq/Mp. Furthermore, there exists a 
parameter region which results the proton life time consistent with the experimental 
results. 
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1 Introduction 



The renormalization group (RG) analysis seems to favor supersymmetric (SUSY) grand 
unified theories (GUTs) over the non-supersymmetric ones. In particular, with the par- 
ticle contents of the minimal supersymmetric standard model (MSSM), the three gauge 
coupling constants converge at the GUT scale Mq — 2 x 10^^ GeV PIE]- In addition, the 
recent progress in neutrino physics [HI makes S'O(IO) GUTs |3] the favorite candidate for 
grand unified theories because it naturally incorporates the see-saw mechanism that can 
naturally explain the lightness of the light neutrino masses. 

Recently, there has been a lot of attention paid to propose and investigate a "minimal" 
5*0 (10) model with so few super multiplets that it can not only fit the current sets of 
Standard Model data and can even predict a few neutrino-related parameters that exper- 
iments are yet to measured accurately. One example of such minimal S'O(IO) model uses 
the irreducible representations 10 + 126 + 126 in additional to the usual quark and lepton 
multiplets of three 16j {i = 1,2,3) and only renormalizable operators 13 El El • 

One of the main undesirable feature of this approach is that, with the sizes of the 
super multiplets employed, they contribute such high beta function to the RG evolution 
such that the GUT gauge coupling constant very quickly blows up to infinity soon after 
the unification scale, Mq- For example, in the model with a set of Higgs representations 
{10 + 126 + 126 + 210} [ini CH, the coupling constant diverges at 4.2 x Mq. While 
this cannot a priori rule out the model, however it does indicate some unknown physics 
may take over even before we reach the string scale or the Planck scale. One possibility 
to explain this run away coupling constant phenomena is to argue that the string scale is 
actually very near the GUT scale such as in some M-theories ^2]- However, this would 
make the success of GUT-related phenomenology more dubious, since we can not neglect 
non- renormalizable operators originated from unknown new physics just above the GUT 
scale. It may be desirable, if achievable, to keep the GUT coupling constant perturbative 
for at least a couple of order of magnitude before it reaches the (perturbative) string scale 
~ 5 X 10^^ GeV [IS] or the (reduced) Planck scale Mp ~ 2.4 x 10^^ GeV. 

On the other hand, the desert scenario associated with the success of the MSSM cou- 
pling constant unification dictates that the GUT scale has to be only about two orders of 
magnitude lower than the Planck scale. It is unavoidable that some higher dimensional 
operators induced by the higher string scale, Mgt or the Planck scale Mp, may play some 
crucial phenomenological role in the analysis of the GUT models. This of course can make 
the simple GUT models much less predictive. However, it is unnatural also to analyze 
GUT models pretending that the string or the Planck scales are not out there not too 
far away. One reasonable strategy to pursue predictability is to use only a minimal set 
of higher dimensional operators as dictated by the requirement of fitting the low energy 
phenomenology. This bottom up approach will leave it to the eventual string or Planck 
scale physics to explain why only these subsets of higher dimensional operators should play 
important role in the GUT model analysis. 

With these perspective in mind, in this paper, we propose a different approach to the 
5*0(10) unification. We pose the question: is it possible to have realistic 5*0(10) unification 
with perturbative coupling constant up to the Planck scale (or the string scale)? We require 
the GUT model to have: 

(1) The coupling constant unification similar to that of MSSM. This will require that even 
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Table 1: List of the Dynkin index for the 5'O(10) irreducible representations up to the 210 
dimensional one 



if there is intermediate scales below the GUT scale, it will have to quite close to Mq. 

(2) The GUT coupling constant remain perturbativc, say ao < 1, up to the Planck scale. 
While it is possible the new physics may come in at the string scale lower than the Planck 
scale, here, in the first analysis, we use Planck scale because it gives stronger constraint on 
the beta function of the GUT theory. 

(3) The GUT model should fit all known low energy experimental data for the Standard 

Model parameters including CP violating phase. There is an issue of the role played by 
the yet undetermined soft SUSY breaking terms. Here we shall assume initially that they 
play no role in the fit to low energy parameters. It is partly because soft SUSY breaking 
sector is the most uncertain part of this analysis. It is reasonable to leave them out until 
it is determined later that they are needed to perfect the model. 



2 Perturbative 50(10) 

The requirement that the 5*0(10) gauge coupling constant remains perturbative up to Mp 
imposes severe constraint on the set of matter and Higgs representations we can use. To 
derive this constraint, note the solution of the (one-loop) RG equation for the unified gauge 
coupling q;g. 



where b = —hgauge + bmatter + ^Higgs IS the beta function coefficient. Each chiral super 
multiplet contributes 1/2 to b, and each vector (gauge) multiplet contributes 3//2 where / 
is the Dynkin index of the irreducible representation listed in Table 1. For 5*0(10) with 
three families, bgauge = 24 and 

bmatter — 2x3, therefore b — — 18 -|- buiggg- If we take the 
constraint and allows the coupling constant to blow up at /i — A, namely 1/q;g(A) = 0, we 
obtain 

In MSSM RG analysis, one typically finds 1/q;(Mg) ^ 24. Therefore bniggs < 49, if one 
uses A = Mp ~ 2.4 x 10^^ GeV, the reduced Planck mass. If one use the stricter condition 
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that a{Mp) = 1, then the constraint becomes bniggs < 48 which is about the same as 
before. Clearly to keep the couplings perturbative, it is necessary to largely reduces the 
Higgs representations. It is clear from Table 1 that the Higgs representations {126 + 126} 
or 210 are forbidden to be introduced into a model. Here note that 126 is necessary 
if 126 is used to break symmetry because of D-flatness condition needed for preserving 
supersymmetry at the GUT scale. 

3 Classifying models 

Two main tasks of the Higgs representations are (1) to break the 5*0(10) gauge symmetry 
down to the Standard Model one and (2) to give fermion masses and mixings consistent 
with all the current experimental data. While there are a priori many choices for Higgs 
representations, we may pick up some Higgs representations to make our model as simple 
as possible. 

For gauge symmetry breaking, the possible choices are: 
(a) {45+54} which contribute bmggs = 20; (b) {16+16+45} which contributes bniggs = 12; 
(c) {16 + 16 + 54} which contributes bniggs — 16. 

The possibility (a) has been analyzed in the literature before. It certainly achieve the 
symmetry breaking down to Gsm — SU (3) x SU (2) x . For (b) , a simple superpotential, 

W = M^545h^5h + Mi6T6ijl6if + XT6h16h45h (3) 

achieves the first task, and one can show that in this case, 45h can develop VEV that 
breaks SO{10) G2231 = SU{2) x SU{2) x SU{3) x U{1) and Hie can further break it to 
Gsm- Note that, in this case, since the usual quark and lepton multiplets also belong to 
a 16, it is necessary to impose a global symmetry like i?-parity to distinguish between the 
usual matter and Higgs multiplets. For (c), the superpotential is only 

W = M5^54h54h + Mi6T6h16h . (4) 

Therefore symmetry breaking is not possible. However, there are still the potential of using 
the higher dimensional operators to help with symmetry breaking. We shall not treat this 
more complicated possibility in this manuscript. 

For fermions masses, there are lots of possible choices for Higgs representations and the 
situation can be more complex. We will define our "minimal model" as the one that can 
accomplish all the above tasks and contributes bniggs as small as possible. 

Considering that top Yukawa coupling is of order one, it is necessary to introduce, 
at least, one Higgs representation which has a renormalizablc Yukawa coupling with 16 
matters. Although both of 10 and 126 can accomplish this task, 126 Higgs is forbidden 
as discussed in the previous section. Thus, we introduce one 10 Higgs into our model. 
Moreover, in order to incorporate Majorana masses of right-handed neutrinos, 16 Higgs is 
necessary and the superpotential 

r|l6,16,T6^T6^ , (5) 

can provide Majorana masses of right-handed neutrinos through VEVs of the MSSM singlet 
components in 16 h and 16 h- It leads to a natural scale of the right-handed neutrino mass 
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as the one derived from the neutrino oscillation data, 

Mr^^. (6) 

Mp ^ ' 

Throughout this paper, we assume only these MSSM singlet components of 16// and IQh 
develop their VEVs. This assumption is essential to write down the GUT mass matrix 
relations for the charged fermions (see the next section). Then, the most reasonable choice 
of the Higgs representations would be {10 + 16 + 16 + 45}. 

With these Higgs representations, the superpotential possibly relevant to the fermion 
masses is given by (up to dimension 5 terms) 

W = r/^16a6,10H + ^n^6a6,10H45// + ^F|l646,T6,/T6,,, (7) 

where the Yukawa coupling matrices Yio, are symmetric, while 145 is antisymmetric. 
Here we have omitted a term proportional to 16jl6jl6/fl6//, since this is irrelevant to the 
fermion mass matrix under the above assumption. One can introduce some global symmetry 
(Zjv symmetry, for example) to forbid some superpotential terms from the beginning, so 
that the model becomes simpler. Such global symmetry also plays a crucial role to forbid 
some dimension five operators (such as 16jl6jl6;;:16i/Mp) in the starting Lagrangian, 
which causes too rapid proton decay. 

In the second term in Eq. ((7j), a product 10//45/f plays the same role as a 120 Higgs 
representation. After VEVs of the Higgs doublets in IQh and 45// in the B — L direction 
are developed, the first two terms in Eq. ((7j) provide Dirac mass matrices of quarks and 
leptons. Note however that this model so far is obviously unrealistic since it predicts the 
Kobayashi-Maskawa matrix being unity. This is because the two terms can be factorized 
by the same IQh and, as result, the up-type quark mass matrix is proportional to the 
down-type quark mass matrix. A simple way to ameliorate this problem is to introduce a 
new 10 Higgs and the superpotential such as 

W = Fi^l6a6,10i + YlllQdQjlQ2^^H , (8) 

where lOi and IO2 are two Higgs multiplets of 10 representation. Here one can again 
introduce some global symmetry under which lOi, IO2 and 45 transform differently, so that 
the couplings of lOi and IO2 are arranged as above. In this case, the second term plays the 
same role of the elementary Higgs of 120 representation and thus this system is effectively 
the same as the one with 10 + 120 elementary Higgs multiplets. Then, our "minimal 
model" is defined by the choice of the set of Higgs representations {2 x 10+16 + 16 + 45}^. 

4 Fermion mass matrices and low energy data fitting 

In the following, we use effective 120 Higgs representation in the analysis. The Yukawa 
couplings relevant to the Dirac mass matrices are given by 

W = YiilQdQjlQn + F/^ol6il6,120// , (9) 



^The "minimal" models similar to our model has been proposed by many authors |14j . 
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where Yiq and Yuo are symmetric and anti-symmetric, respectively. Here note that Y120 = 
145 (45//) /Mp in the terms of the original superpotential with the VEV of 45 h in the 
B — L direction. Both of the Higgs multiplets IOh and 120// include a pair of Higgs 
doublets in the MSSM decomposition. At low energy after the GUT symmetry breaking, 
the superpotential leads to 

+ {Yl^Hl, - ?,Yli,Hl^,)Nd, + {YllHt, - ?,Yli,HU)e% , (10) 

where if 10 and H120 correspond to the Higgs doublets in 10// an 120//, which originate lOi 
and IO2 in the superpotential of Eq. (jSJ. The factor —3 in the lepton sector is the results 
from the VEV of 45// in the B ~ L direction, and plays a crucial role so that unwanted 
GUT mass relations, me = irtd and = rris, is corrected. 

In order to keep the successful gauge coupling unification, suppose that one pair of 
Higgs doublets given by a linear combination of Hiq'^ and H120 is light while the other pair 
is heavy (~ Mg)- The light Higgs doublets are identified as the MSSM Higgs doublets {Hu 
and Hd) and given by 

Hd = OidH^Q + (3dHi2o , (11) 

where au,d and jSu^d denote elements of the unitary matrix which rotate the flavor basis 
in the original model into the (SUSY) mass eigenstates. Omitting the heavy Higgs mass 
eigenstates, the low energy superpotential is described by only the light Higgs doublets Hu 
and Hd such that 

Wy = u^{a^Yll+P^Yli,)Huq^ + d^\a''Yii + P''Yii,)Hdq, 

+ (a^Yli - WYlQ Hu + e\ (a'^YH - SP'^Y^Q H^ i, , (12) 

where the formulas of the inverse unitary transformation of Eq. (Ill|) . Hi^f" = a^^'^H^^d + ■ ■ ■ 
and H120 = P^''^Hu4 + ■ ■ ■, have been used. 

Providing the Higgs VEVs, {Hy) = f sin/? and {Hd) = vcosP with v ~ 174 GeV, the 
Dirac mass matrices can be read off as 

Mu = cioMio + C120M120 , 
Md = M10 + M120, 

Md = CioMio - 3Ci2oMl20 , 

Me = M10-3M120, (13) 

where M„, Md, Md and Me denote up-type quark, down- type quark, neutrino Dirac, 
charged-lepton mass matrices, respectively. Note that all the mass matrices are described 
by using only two basic mass matrices, a symmetric Mio and an antisymmetric M120, and 
two complex coefficients Cio and C120, which are defined as Miq = Yi^a'^vcosP, M120 = 
Yi2o/3'^vcos/3, cio = {a^ /a^) tan l3 and C120 = {/3'^ //3'^) tan l3, respectively. 

These mass matrix formulas lead to the GUT mass matrix relation among the quark 
and lepton mass matrices, 

M, = CdiMd + KMu) , (14) 
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where 



3cio + C120 
Cd = , 

ClO — C120 

4 

n = . (15) 

3cio + C120 

For simplicity, we assume that Mio and M120 are real and pure imaginary matrices, re- 
spectively, and Cio and C120 are both real. Then, all the Dirac mass matrices becomes 
hermitian and still CP violating. Note that, according to this assumption, the num- 
ber of free parameters in our model are reduced into eleven in total; six real parameters 
in Mio, three real parameters in M120, Cio and Ci2o- On the other hand, the number of 
observables we should fit is thirteen; six quark masses, three angles and one CP-phase in 
the CKM matrix and three charged lepton masses. Thus there are two predictions for 
observables, whose values have been already known by experiments. Therefore, the data 
fitting in our model is very non-trivial. In the following analysis, the strange quark mass 
and the CP-phase in the CKM matrix will be regarded as two predictions in our model 
(see the following discussion). 

Without loss of generality, we can begin with the basis where is real and diagonal, 
Mu = Du- In this basis, the hermitian matrix can be described as = VcKhtDdVcKM 
by using the CKM matrix Vckm and the real diagonal mass matrix ■ Considering the 
basis-independent quantities, tr (Me), tr (M^) and det (Mg), and eliminating q, we obtain 
two independent equations. 



tr(Me) V tr(M/) 



me + -|- rrir ) mj + -|- 

tr{Me) \^ det (Me) 



(16) 
(17) 



\me + + rrir J me 

where Me = Vckm Dd V^t^^j -|- kDu- With input data of six quark masses, three angles and 
one CP-phase in the CKM matrix and three charged lepton masses, we solve the above 
equations and determine k. Using the k determined, Cd is fixed by 

^ _ me + + TTlr 

{rrid + nis + m^) -|- k {m^ + m^ + m^) 
The original basic mass matrices, Mio and M120, are described by 

Mio = ^^VcKMDdV^KM + ''-^D^. (19) 
M120 = ^^VcKwiDdV^KM-^^D^. (20) 

Once the solutions q and k are obtained, Mio and M120 are completely determined. 

Note that it is a very non-trivial problem to find a solution that satisfies both Eqs. (fT^ 
and ()17|1 at the same time with only one free parameter k. In the following analysis, we vary 

^In general, Md — U DdU^ by using a general unitary matrix U = e*"e*'^"^^e*'''^** Vci^AfC*'' '^^e^^ In 
this paper, we adopt the diagonal phases to or tt for simplicity. 
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two input data, the strange quark mass (mg) and the CP-phase [S) in the CKM matrix, 
within their experimental errors, so that both Eqs. ()16|) and fjl7|) can be satisfied with the 
same k, value in good accuracy. We can find a consistent solution only if we input special 
values for and 6. This fact indicates that the input values for and S we have used 
in our analysis are two predictions in our model as mentioned above. 

Now let us solve the GUT relation and determine Cd and k. We follow the same strategy 
in [Hj. Since the GUT mass matrix relation is valid only at the GUT scale, we first evolve 
the data at the weak scale to the ones at the GUT scale with given tan (3 according to the 
renormalization group equations (RGEs) and use them as input data at the GUT scale. 
We take input the absolute values of the fermion masses at Mz as follows (in GeV): 

ruu = 0.00233, m, = 0.677, = 176, 
ma = 0.00469, = 0.0747, nib = 3.00, 

me = 0.000487, = 0.103, = 1.76. (21) 

Here the experimental values extrapolated from low energies to Mz were used ^3], and we 
choose the signs of the input fermion masses as (m„, rric, rrit) = (— , — , +) and (m^, m^, rrih) = 
(—,—,+). For the CKM mixing angles and a CP- violating phase in the "standard" param- 
eterization, we input the values measured by experiments as follows: 

si2 = 0.2229, S23 = 0.0412, = 0.0036, 5 = 57.6° . (22) 

Since it is very difficult to search all the possible parameter region systematically, we present 
our results for tan/3 = 30. Note that only the case of large tan/3 can be consistent with 
our original model, since 1^20 = ^45(45)/Mp ~ 0.011^*5 and only y/J can be of order one. 
After the RGE running, we obtain the fermion masses and the CKM mixing angles and the 
CP phase at the GUT scale, and use them as input parameters in order to solve Eqs. (fTB|) 
and (fTTj) . By putting the above data, Eq. (fT^ gives a solution 

K = -0.01107011- ■■ . (23) 

On the other hand, from Eq. p7|) . we obtain the solution 

K = -0.01107006- ■ ■ . (24) 

Since these solutions are coincide with each other in good accuracy, we can regard it as the 
solution we seek. Using Eq. (fT^ to determine q, now we find a solution 

K = -0.0111 , 

Cd = -7.89 . (25) 

The existence of the solution means that our model can reproduce the low energy experi- 
mental data for the charged fermion sector. 

5 Proton decay 

The most characteristic prediction of the SUSY GUTs is the proton decay. Normally 
in SUSY GUTs the proton decay process through the dimension five operators involving 
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MSSM matters, mediated by the color triplet Higgsino, turns out to be the dominant 
decay modes, since the process is suppressed by only a power of the Higgsino mass scale. 
Experimental lower bound on the proton decay modes p — *• K^V through the dimension 
five operators is given by SuperKamiokande (SuperK) P^ . 

rip K+V) > 2.2 X 10^^ [years] . (26) 

This is one of the most stringent constraints in construction of phenomenologically viable 
SUSY GUT models. In fact, the minimal SUSY SU{5) model has been argued to be 
excluded from the experimental bound together with the requirement of the success of the 
three gauge coupling unification I17 j . However note that the minimal SU{5) model 
predictions contradicts against the realistic charged fermion mass spectrum, and thus, 
strictly speaking, the model is ruled out from the beginning. Obviously some extensions of 
the fiavor structures in the model is necessary to accommodate the realistic fermion mass 
spectrum. On the other hand, knowledge of the fiavor structure is essential in order to give 
definite predictions about the proton decay processes through the dimension five operators. 
Some models in which fiavor structures are extended have been found to be consistent with 
the experiments in the context of SU{5) models [T^ 1^ and in 5*0(10) extensions ^T] . 

As discussed in the previous section, the charged fermion mass matrices have been 
completely determined in our model. Therefore, we can investigate proton decay rate with 
only some free parameters. Our discussion follows [221 • 

The Yukawa interactions of the MSSM matter with the color triplet Higgs induces the 
following Baryon and Lepton number violating dimension five operator 

W = C'i^^Q'Q^Q^L^ . (27) 

Here the coefficients are given by the products of the Yukawa coupling matrices and the 
(effective) color triplet Higgsino mass matrix. In our model, the coefficients are given by 
the products of two basic Yukawa coupling matrices, Yio and Y120, and the effective 2x2 
color triplet Higgsino mass matrix, Mc, such as [TT] 



Cr = (Yli, n^o) (Mc') ( ^fi ) . (2^ 

120 / 



Vkl 

yk 
1 1 ' 



As discussed in the previous section, the Yukawa coupling matrices, Fio and F120, are 
related to the corresponding mass matrices Mio and M120 such as 



a^v sm p 



Yuo = -^^M,2o, (29) 
p"f sm p 

with f ~ 174 GeV. Here a" and are the Higgs doublet mixing parameters introduced 
in the previous section, which are restricted in the range la**!^ + |/?"|^ < 1. Although these 
parameters are irrelevant to fit the low energy experimental data of the charged fermion 
mass matrices, there is a theoretical lower bound on them in order for the resultant Yukawa 
coupling constant not to exceed the perturbative regime. In order to obtain the most con- 
servative values of the proton decay rate, we make a choice of the Yukawa coupling matrices 



8 



as small as possible. In the following analysis, we restrict the region of the parameters in 
the range [a^)"^ + = 1 (we assume a" and /3" real for simphcity). Here we present an 
example of the Yukawa coupling matrices with fixed a" = 0.202 (/5" = 0.979), 

/ 0.00261 0.00485 -0.00208 \ 
^10 = 0.00485 0.0163 -0.0414 , (30) 
V -0.00208 -0.0414 1.00 / 



120 



-0.0000379 i -0.0000379 i 

0.0000379 z -4.61 X 10-S | . (31) 

V 0.0000379 z 4.61 x 10~^ i 



Note that the numerical smallness of Y120 is a consistency check of our scheme. Since it is 
a result of the higher dimensional operator, its smallness relative to indicated that it is 
reasonable to ignore the even higher dimensional operator. However, Y'120 itself does play 
important role in fitting low energy data including CP violation. 

For the effective color triplet Higgsino mass matrix, we assume degenerate eigenvalues 
being the GUT scale, Mq = 2 x 10^^ GeV, which is necessary to keep the successful gauge 
coupling unification. Then, in general, we can parameterize the 2x2 mass matrix as 

Mc = MgxU , (32) 

with the unitary matrix, 

U^^i^^rJ cose sme ^e'^'"^. (33) 
\ — sm u cos I 

Here we omit an over all phase since it is irrelevant to calculations of the proton decay rate. 
There are four free parameters in total involved in the coefficient C^j^''\ namely, a", ip, ip' 
and 6. Once these parameters are fixed, (7^^^' is completely determined. 

Through the same numerical analysis as in [22j we can find a parameter region in which 
the proton life time can be in the range consistent with the SuperK results. In fact, we can 
find a special colored Higgsino mass matrix that can cancel the proton decay rate through 
the dominant mode p K^Vt-. For example, for the Yukawa coupling matrices of Eqs. (jHUj) 
and (jSD), it is found to be (tan/3 = 30) 

( -0.0681i 0.998 \ 

in other words, 6 = 1.64 [rad], ip = 1.57 [rad], and ip' = [rad]. With these parameters 
fixed, the proton life time through the sub-dominant decay modes is estimated as follows: 

r{p K^Ve) = 2.5 X 10^^ [years] , (35) 
t{p K^V^) = 4.3 X 10=^=^ [years] . (36) 

In our analysis, we have taken the averaged squark mass of the 1st and 2nd generations 
as m~ = 10 TeV and the Wino mass as = 500 GeV. These results exceed the current 
experimental lower bounds. Therefore, our model passes the proton decay constraint and 
is phenomenologically viable. 
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6 Neutrino physics 



In our model, the right-handed Majorana neutrino mass matrix is generated by Eq. (0) 
through VEV of 16^ in the MSSM singlet direction. Here Yjg is the complex symmetric 
matrix which has twelve free parameters in general and it is nothing to do with charged 
fermion mass matrices. Therefore, through the see-saw mechanism, there is enough number 
of free parameters to fit all the current neutrino oscillation data. In other words, there is no 
prediction for neutrino oscillation physics. However, there is an interesting feature through 
the see-saw mechanism. 

In the basis where Mg is (positive) real and diagonal, the light neutrino mass matrix 
is given by the see-saw mechanism Mj, = M^M^^M/j. M^, is diagonalized by the Maki- 
Nakagawa-Sakata (MNS) mixing matrix such as M^, = Uj,j^gdia.g{mi,m2,m^)UMNS- The 
current neutrino oscillation data provide informations (but not complete) for the mixing 
angles 9ij in the MNS matrix and m,. 

Recall that, as discussed in Sec. 4, all the elements in the neutrino Dirac mass matrix 
can be determined in our model. Therefore, information of can be extracted through 
the (inverse) see-saw relation. 



if the neutrino oscillation data are used as inputs. Since the current experimental data are 
insufficient to fix all the elements in My, Mji can be described as a function of parameters 
not yet undetermined by experiments, Mr = Mij(m£, ^13, 5, /?, 7), where me is the lightest 
mass eigenvalues of the light Majorana neutrino, 6, (3 and 7 are the Dirac CP-phase and the 
Majorana CP-phases, respectively. Making some assumptions for these free parameters, 
one can evaluate Mr concretely, and leads to predictions for physics related to the right- 
handed neutrino mass matrix, such as, the leptogenesis scenario [22] • This direction would 
be worth investigating. We leave it for future works. 

In addition, an order estimation leads to an interesting consequence. In our model, the 
neutrino Dirac mass matrix is approximately the same as the up-type quark mass matrix, 
and hence its heaviest eigenvalue is roughly the same as the top quark mass. As already 
mentioned, the natural scale of the right-handed neutrino mass is M^ ~ Mq/Mp ~ 10^^ 
GeV. Therefore, according to the see-saw mechanism, we find the heaviest mass eigenvalue 



of the light neutrinos being of order 0.1 eV. Interestingly, this value is close to yAm^, 

where Am^ ~ 2.1 x 10"'^ eV^ is the atmospheric neutrino oscillation data [3]. This result 
indicates that our model prefers the hierarchical case to the degenerate case for the light 
neutrino mass spectrum. 

7 Conclusion 

Discovery of the neutrino masses and mixings has made 50(10) GUT models the favorite 
candidate as new physics. Lots of 5*0(10) GUT models have been intensively discussed. 
There are a priori many choices for the Higgs representations to be introduced into a 
model. We have imposed the requirement that the GUT gauge coupling should remain 
to be perturbative up to the (reduced) Planck scale or the string scale at which further 
new physics including quantum gravity takes over. This requirement has been found to 



Mr 



Md M-' Ml , 



(37) 




10 



be strong enough to forbids Higgs representations higher than 126 dimension. We have 
proposed a simple 5*0(10) model with a set of Higgs representations {2x 10 + 16+16 + 45}, 
which can satisfy the requirement and gives the beta function coefficient of the GUT gauge 
couphng as small as possible. It has been shown that the model is phenomenologically 
viable, namely all the charged fermion masses and mixings have been reproduced. In this 
realistic Yukawa couplings, the most stringent proton decay processes has been suppressed. 
Also, the model can reproduce the current neutrino oscillation data with the right-handed 
neutrino mass being of order Mr ~ Mq/Mp. 
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